We study data on Mongolia and Japan.

1 WDI

knitr::include_url("https://icu-hsuzuki.github.io/science/index-j.html")
$url
[1] "https://icu-hsuzuki.github.io/science/index-j.html"

$height
[1] "400px"

$url.orig
[1] "https://icu-hsuzuki.github.io/science/index-j.html"

attr(,"class")
[1] "knit_embed_url" "knit_asis_url" 

2 Setup

Sys.setenv(LANG = "en")
library(tidyverse)
library(WDI)

3 Education

3.1 Education in Japan

url_edu1 <- "https://www.gender.go.jp/about_danjo/whitepaper/r03/zentai/html/honpen/csv/zuhyo01-05-01.csv"
download.file(url_edu1, destfile = "data/edu1.csv")
trying URL 'https://www.gender.go.jp/about_danjo/whitepaper/r03/zentai/html/honpen/csv/zuhyo01-05-01.csv'
Content type 'text/csv' length 3220 bytes
==================================================
downloaded 3220 bytes
guess_encoding("data/edu1.csv")
df_edu1 <- read_csv("data/edu1.csv", locale = locale(encoding = "Shift-JIS"), skip = 2)
Rows: 71 Columns: 10
── Column specification ───────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (1): 年度
dbl (9): 高等学校等(男子), 高等学校等(女子), 専修学校(専門課程,男子), 専修学校(専門課程,女子), 大学(学部,男子), 大学(学部,女子), 短期大学(本科,女子), 大学院...

ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.
df_edu1
colnames(df_edu1)
 [1] "年度"                       "高等学校等(男子)"         "高等学校等(女子)"        
 [4] "専修学校(専門課程,男子)" "専修学校(専門課程,女子)" "大学(学部,男子)"        
 [7] "大学(学部,女子)"         "短期大学(本科,女子)"     "大学院(男子)"            
[10] "大学院(女子)"            
df_edu1$年度
 [1] "昭和25" "26"     "27"     "28"     "29"     "30"     "31"     "32"     "33"     "34"     "35"    
[12] "36"     "37"     "38"     "39"     "40"     "41"     "42"     "43"     "44"     "45"     "46"    
[23] "47"     "48"     "49"     "50"     "51"     "52"     "53"     "54"     "55"     "56"     "57"    
[34] "58"     "59"     "60"     "61"     "62"     "63"     "平成元" "2"      "3"      "4"      "5"     
[45] "6"      "7"      "8"      "9"      "10"     "11"     "12"     "13"     "14"     "15"     "16"    
[56] "17"     "18"     "19"     "20"     "21"     "22"     "23"     "24"     "25"     "26"     "27"    
[67] "28"     "29"     "30"     "令和元" "2"     
df_edu0 <- df_edu1
colnames(df_edu0) <- c("year", "highschool_m", "highschool_f", "vocational_m", "vocational_f", "university_m", "university_f", "juniorcol_f", "gradschool_m", "gradschool_f")
df_edu00 <- df_edu0 %>% mutate(year = 1950:2020, 
                   highschool = (highschool_m + highschool_f)/2,
                   vocational = (vocational_m + vocational_f)/2,
                   university = (university_m + university_f)/2, 
                   juniorcol = juniorcol_f,
                   gradschool = (gradschool_m + gradschool_f)/2) 
df_edu00 %>% select(-(2:10)) %>%
  pivot_longer(!year, names_to = "schools", values_to = "percentage") %>%
  ggplot(aes(x = year, y = percentage, color = schools)) +
  geom_line()
Warning: Removed 49 row(s) containing missing values (geom_path).

df_edu11 <- df_edu1 %>% mutate(年 = 1950:2020, 
                   高等学校 = (`高等学校等(男子)`+ `高等学校等(女子)`)/2,
                   専修学校 = (`専修学校(専門課程,男子)`+ `専修学校(専門課程,女子)`)/2,
                   大学 = (`大学(学部,男子)` + `大学(学部,女子)`)/2, 
                   短期大学 = `短期大学(本科,女子)`,
                   大学院 = (`大学院(男子)` + `大学院(女子)`)/2) %>%
  select(-(1:10))
df_edu11
df_edu11 %>% 
  pivot_longer(2:6, names_to = "学校", values_to = "進学率") %>%
  ggplot(aes(x = 年, y = 進学率, color = 学校)) +
  geom_line()
df_edu11 %>% 
  pivot_longer(2:6, names_to = "学校", values_to = "進学率") %>%
  ggplot(aes(x = 年, y = 進学率, color = 学校)) +
  geom_line() +
  theme_gray (base_family = "HiraginoSans-W3") # or base_family = "HiraKakuPro-W3"
Warning: Removed 49 row(s) containing missing values (geom_path).

df_edu1 %>% mutate(year = 1950:2020) %>%
  pivot_longer(2:10, names_to = "学校", values_to = "進学率") %>%
  ggplot(aes(x = year, y = 進学率, color = 学校)) +
  geom_line() +
  theme_gray (base_family = "HiraKakuPro-W3") # or base_family = "HiraginoSans-W3"
Warning: Removed 94 row(s) containing missing values (geom_path).

#{r dev='rag_png'} df_edu1 %>% mutate(year = 1950:2020) %>% pivot_longer(2:10, names_to = "学校", values_to = "進学率") %>% ggplot(aes(x = year, y = 進学率, color = 学校)) + geom_line()

LS0tCnRpdGxlOiAiTW9uZ29saWEgYW5kIEphcGFuIgphdXRob3I6ICJEUy1TTCIKZGF0ZTogImByIFN5cy5EYXRlKClgIgpvdXRwdXQ6CiAgaHRtbF9ub3RlYm9vazogCiAgICBudW1iZXJfc2VjdGlvbnM6IHllcwogICAgdG9jOiB0cnVlCiAgICB0b2NfZmxvYXQ6IHRydWUKICBodG1sX2RvY3VtZW50OgogICAgZGZfcHJpbnQ6IHBhZ2VkCi0tLQo+IFdlIHN0dWR5IGRhdGEgb24gTW9uZ29saWEgYW5kIEphcGFuLgoKIyBXREkKCiogaHR0cHM6Ly9kYXRhLndvcmxkYmFuay5vcmcvP2xvY2F0aW9ucz1NTi1KUC0xVwogIC0gQ29tcGFyaXNvbiBvZiBkYXRhIG9mIE1vbmdvbGlhLCBKYXBhbiBhbmQgdGhlIFdvcmxkCiAgCjxpZnJhbWUgc3JjPSJodHRwczovL291cndvcmxkaW5kYXRhLm9yZy9ncmFwaGVyL2NoaWxkcmVuLXBlci13b21hbi1VTiIgbG9hZGluZz0ibGF6eSIgc3R5bGU9IndpZHRoOiAxMDAlOyBoZWlnaHQ6IDYwMHB4OyBib3JkZXI6IDBweCBub25lOyI+PC9pZnJhbWU+Cgo8aWZyYW1lIHNyYz0iaHR0cHM6Ly9kYXRhLndvcmxkYmFuay5vcmcvc2hhcmUvd2lkZ2V0P2luZGljYXRvcnM9U1AuUE9QLkdST1cmbG9jYXRpb25zPU1OLUpQLTFXIiBsb2FkaW5nPSJsYXp5IiBzdHlsZT0id2lkdGg6IDEwMCU7IGhlaWdodDogNjAwcHg7IGJvcmRlcjogMHB4IG5vbmU7Ij48L2lmcmFtZT4KICAKPGlmcmFtZSBzcmM9Imh0dHBzOi8vZGF0YS53b3JsZGJhbmsub3JnL3NoYXJlL3dpZGdldD9pbmRpY2F0b3JzPVNQLlBPUC5HUk9XJmxvY2F0aW9ucz1NTi1KUC0xVyIgd2lkdGg9JzMwMCcgaGVpZ2h0PScyMDAnIGZyYW1lQm9yZGVyPScwJyBzY3JvbGxpbmc9Im5vIiA+PC9pZnJhbWU+CiAgCmBgYHs9aHRtbH0KPGlmcmFtZSBzcmM9Imh0dHBzOi8vZGF0YS53b3JsZGJhbmsub3JnL3NoYXJlL3dpZGdldD9pbmRpY2F0b3JzPVNQLkRZTi5MRTAwLklOJmxvY2F0aW9ucz1NTi1KUC0xVyIgd2lkdGg9JzQ1MCcgaGVpZ2h0PSczMDAnIGZyYW1lQm9yZGVyPScwJyBzY3JvbGxpbmc9Im5vIiA+PC9pZnJhbWU+CmBgYAoKYGBge3J9CmtuaXRyOjppbmNsdWRlX3VybCgiaHR0cHM6Ly9pY3UtaHN1enVraS5naXRodWIuaW8vc2NpZW5jZS9pbmRleC1qLmh0bWwiKQpgYGAKCiMgU2V0dXAKCmBgYHtyIHNldHVwfQpTeXMuc2V0ZW52KExBTkcgPSAiZW4iKQpsaWJyYXJ5KHRpZHl2ZXJzZSkKbGlicmFyeShXREkpCmBgYAoKIyBFZHVjYXRpb24KCiMjIEVkdWNhdGlvbiBpbiBKYXBhbgoKYGBge3IgY2FzaCA9IFRSVUV9CnVybF9lZHUxIDwtICJodHRwczovL3d3dy5nZW5kZXIuZ28uanAvYWJvdXRfZGFuam8vd2hpdGVwYXBlci9yMDMvemVudGFpL2h0bWwvaG9ucGVuL2Nzdi96dWh5bzAxLTA1LTAxLmNzdiIKZG93bmxvYWQuZmlsZSh1cmxfZWR1MSwgZGVzdGZpbGUgPSAiZGF0YS9lZHUxLmNzdiIpCmd1ZXNzX2VuY29kaW5nKCJkYXRhL2VkdTEuY3N2IikKYGBgCgpgYGB7cn0KZGZfZWR1MSA8LSByZWFkX2NzdigiZGF0YS9lZHUxLmNzdiIsIGxvY2FsZSA9IGxvY2FsZShlbmNvZGluZyA9ICJTaGlmdC1KSVMiKSwgc2tpcCA9IDIpCmRmX2VkdTEKYGBgCgpgYGB7cn0KY29sbmFtZXMoZGZfZWR1MSkKZGZfZWR1MSTlubTluqYKYGBgCgpgYGB7cn0KZGZfZWR1MCA8LSBkZl9lZHUxCmNvbG5hbWVzKGRmX2VkdTApIDwtIGMoInllYXIiLCAiaGlnaHNjaG9vbF9tIiwgImhpZ2hzY2hvb2xfZiIsICJ2b2NhdGlvbmFsX20iLCAidm9jYXRpb25hbF9mIiwgInVuaXZlcnNpdHlfbSIsICJ1bml2ZXJzaXR5X2YiLCAianVuaW9yY29sX2YiLCAiZ3JhZHNjaG9vbF9tIiwgImdyYWRzY2hvb2xfZiIpCmRmX2VkdTAwIDwtIGRmX2VkdTAgJT4lIG11dGF0ZSh5ZWFyID0gMTk1MDoyMDIwLCAKICAgICAgICAgICAgICAgICAgIGhpZ2hzY2hvb2wgPSAoaGlnaHNjaG9vbF9tICsgaGlnaHNjaG9vbF9mKS8yLAogICAgICAgICAgICAgICAgICAgdm9jYXRpb25hbCA9ICh2b2NhdGlvbmFsX20gKyB2b2NhdGlvbmFsX2YpLzIsCiAgICAgICAgICAgICAgICAgICB1bml2ZXJzaXR5ID0gKHVuaXZlcnNpdHlfbSArIHVuaXZlcnNpdHlfZikvMiwgCiAgICAgICAgICAgICAgICAgICBqdW5pb3Jjb2wgPSBqdW5pb3Jjb2xfZiwKICAgICAgICAgICAgICAgICAgIGdyYWRzY2hvb2wgPSAoZ3JhZHNjaG9vbF9tICsgZ3JhZHNjaG9vbF9mKS8yKSAKYGBgCgpgYGB7cn0KZGZfZWR1MDAgJT4lIHNlbGVjdCgtKDI6MTApKSAlPiUKICBwaXZvdF9sb25nZXIoIXllYXIsIG5hbWVzX3RvID0gInNjaG9vbHMiLCB2YWx1ZXNfdG8gPSAicGVyY2VudGFnZSIpICU+JQogIGdncGxvdChhZXMoeCA9IHllYXIsIHkgPSBwZXJjZW50YWdlLCBjb2xvciA9IHNjaG9vbHMpKSArCiAgZ2VvbV9saW5lKCkKYGBgCmBgYHtyfQpkZl9lZHUxMSA8LSBkZl9lZHUxICU+JSBtdXRhdGUo5bm0ID0gMTk1MDoyMDIwLCAKICAgICAgICAgICAgICAgICAgIOmrmOetieWtpuagoSA9IChg6auY562J5a2m5qCh562J77yI55S35a2Q77yJYCsgYOmrmOetieWtpuagoeetie+8iOWls+WtkO+8iWApLzIsCiAgICAgICAgICAgICAgICAgICDlsILkv67lrabmoKEgPSAoYOWwguS/ruWtpuagoe+8iOWwgumWgOiqsueoi++8jOeUt+WtkO+8iWArIGDlsILkv67lrabmoKHvvIjlsILploDoqrLnqIvvvIzlpbPlrZDvvIlgKS8yLAogICAgICAgICAgICAgICAgICAg5aSn5a2mID0gKGDlpKflrabvvIjlrabpg6jvvIznlLflrZDvvIlgICsgYOWkp+Wtpu+8iOWtpumDqO+8jOWls+WtkO+8iWApLzIsIAogICAgICAgICAgICAgICAgICAg55+t5pyf5aSn5a2mID0gYOefreacn+Wkp+Wtpu+8iOacrOenke+8jOWls+WtkO+8iWAsCiAgICAgICAgICAgICAgICAgICDlpKflrabpmaIgPSAoYOWkp+WtpumZou+8iOeUt+WtkO+8iWAgKyBg5aSn5a2m6Zmi77yI5aWz5a2Q77yJYCkvMikgJT4lCiAgc2VsZWN0KC0oMToxMCkpCmRmX2VkdTExCmBgYApgYGAKZGZfZWR1MTEgJT4lIAogIHBpdm90X2xvbmdlcigyOjYsIG5hbWVzX3RvID0gIuWtpuagoSIsIHZhbHVlc190byA9ICLpgLLlrabnjociKSAlPiUKICBnZ3Bsb3QoYWVzKHggPSDlubQsIHkgPSDpgLLlrabnjocsIGNvbG9yID0g5a2m5qChKSkgKwogIGdlb21fbGluZSgpCmBgYAoKYGBge3J9CmRmX2VkdTExICU+JSAKICBwaXZvdF9sb25nZXIoMjo2LCBuYW1lc190byA9ICLlrabmoKEiLCB2YWx1ZXNfdG8gPSAi6YCy5a2m546HIikgJT4lCiAgZ2dwbG90KGFlcyh4ID0g5bm0LCB5ID0g6YCy5a2m546HLCBjb2xvciA9IOWtpuagoSkpICsKICBnZW9tX2xpbmUoKSArCiAgdGhlbWVfZ3JheSAoYmFzZV9mYW1pbHkgPSAiSGlyYWdpbm9TYW5zLVczIikgIyBvciBiYXNlX2ZhbWlseSA9ICJIaXJhS2FrdVByby1XMyIKYGBgCgoKYGBge3J9CmRmX2VkdTEgJT4lIG11dGF0ZSh5ZWFyID0gMTk1MDoyMDIwKSAlPiUKICBwaXZvdF9sb25nZXIoMjoxMCwgbmFtZXNfdG8gPSAi5a2m5qChIiwgdmFsdWVzX3RvID0gIumAsuWtpueOhyIpICU+JQogIGdncGxvdChhZXMoeCA9IHllYXIsIHkgPSDpgLLlrabnjocsIGNvbG9yID0g5a2m5qChKSkgKwogIGdlb21fbGluZSgpICsKICB0aGVtZV9ncmF5IChiYXNlX2ZhbWlseSA9ICJIaXJhS2FrdVByby1XMyIpICMgb3IgYmFzZV9mYW1pbHkgPSAiSGlyYWdpbm9TYW5zLVczIgpgYGAKYGBgI3tyIGRldj0ncmFnX3BuZyd9CmRmX2VkdTEgJT4lIG11dGF0ZSh5ZWFyID0gMTk1MDoyMDIwKSAlPiUKICBwaXZvdF9sb25nZXIoMjoxMCwgbmFtZXNfdG8gPSAi5a2m5qChIiwgdmFsdWVzX3RvID0gIumAsuWtpueOhyIpICU+JQogIGdncGxvdChhZXMoeCA9IHllYXIsIHkgPSDpgLLlrabnjocsIGNvbG9yID0g5a2m5qChKSkgKwogIGdlb21fbGluZSgpCmBgYAoKYGBgez1odG1sfQoKPGlmcmFtZSBzcmM9Imh0dHBzOi8vZGF0YS53b3JsZGJhbmsub3JnL3NoYXJlL3dpZGdldD9pbmRpY2F0b3JzPVNQLkRZTi5MRTAwLklOJmxvY2F0aW9ucz1NTi1KUC0xVyIgd2lkdGg9JzQ1MCcgaGVpZ2h0PSczMDAnIGZyYW1lQm9yZGVyPScwJyBzY3JvbGxpbmc9Im5vIiA+PC9pZnJhbWU+CmBgYAo=